function e=ERP(THETA,T)
    TR = thetaR(T');
    e =(TR-THETA)./TR;
end